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Abstract 



We study soft-gluon radiation for heavy quark production in charged-current 
Deep Inelastic Scattering processes. We resum large- x contributions to the 
MS coefficient function to next-to-leading logarithmic accuracy and present 
results for charm quark structure functions for several values of the ratio 
m c /Q. The effect of soft resummation turns out to be visible, especially at 
small Q 2 , and the results exhibit very little dependence on factorization and 
renormalization scales. The impact of our calculation on structure function 
measurements at NuTeV and HERA experiments is discussed. 



1 Introduction 



Lepton-hadron Deep Inelastic Scattering (DIS) is one of the most interesting processes to 
investigate Quantum Chromodynamics, from both theoretical and experimental point 
of view. In fact, DIS allows one to measure proton structure functions, and hence 
investigate the internal structure of hadrons. For the sake of performing accurate mea- 
surements of observables in DIS processes, precise QCD calculations are necessary. 

In this paper we study heavy quark production in DIS, in particular charm quark 
production in charged-current (CC) events. Measurements of the charm quark structure 
functions in the CC regime are in fact important in order to probe the density of strange 
quarks and gluons in the proton. 

Charged-current Deep Inelastic Scattering processes v^N — > \iX are currently under 
investigation at the NuTeV experiment at Fermilab, where neutrino and antineutrino 
beams collide on an fixed iron target (see, e.g., Ref. P for the updated results). In par- 
ticular, production of two oppositely-charged muons at NuTeV p. is mainly associated 
with a CC event with a charm quark in the final state. 

The first generation of HI and ZEUS experiments on electron-proton DIS at the 
HERA collider at DESY (HERA I) did detect charged-current events ep -> v e X HUE], 
but it did not have sufficiently-high statistics to reconstruct heavy quarks in the CC 
regime. However, such measurements are foreseen in the current improved-luminosity 
run (HERA II) and at the possible future generation of HERA III experiments [H] , which 
should also be able to investigate the proton structure functions at low Q 2 values. 

Structure functions are given by convolutions of short-distance, partonic coefficient 
functions, C(x,fi 2 F ), usually given in the MS factorization scheme, with long-distance 
parton distribution functions f(x,fi 2 7 ), whose dependence on the factorization scale //p 
is ruled by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equa- 
tions [HJIZj- The next-to- leading order (NLO) MS coefficient functions for charm pro- 
duction in CC DIS present terms which get arbitrarily large once the incoming quark 
energy fraction z in qW* — > c(g) approaches unity. This corresponds to soft-gluon radi- 
ation. In order to reliably predict structure functions at large z, such terms need to be 
resummed to all orders in the strong coupling constant as (threshold resummation). It 
is the purpose of this paper to perform the resummation of the MS coefficient function 
and investigate the phenomenological implications of the resummation on charm quark 
structure functions at NuTeV and HERA. 

Soft resummation for the coefficient function in DIS was first considered in |H1I§1ITU] 
and later implemented in but in the approximation in which all participating quarks 
are massless. Soft-gluon resummation for heavy quark production in DIS processes 
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has been investigated in fSlEj, where the authors have considered neutral current 
interactions. Phenomenological studies at NLO on charm quark production at NuTeV 
and HERA have been performed in and ^5] respectively. 

The plan of our paper is the following. In Section 2 we review the NLO result for 
heavy quark production in CC DIS and discuss the behaviour of the coefficient function 
for soft-gluon radiation. In Section 3 we present results for the soft-resummed MS 
coefficient function. In particular, we shall compare the results on CC DIS with soft 
resummation in top quark decay P^ ITTj and we shall comment on the differences with 
respect to the previous work on light quark production in DIS [HI EH HI] • In Section 
4 we present our predictions on charm quark structure functions for several values of 
Q 2 which are relevant for NuTeV and HERA phenomenology. In Section 5 we shall 
summarize the main results of our work and discuss its possible future extensions. 



2 Coefficient functions and structure functions at 
NLO 

In this section we would like to review the main results on CC DIS at NLO in the strong 
coupling constant as. The Born parton-level process in CC DIS reads: 

qi ( Pl )W*(q) - q 2 (p 2 ). (1) 

Process (0) receives 0(as) corrections from quark-scattering 

qi(px)W*(q)^q 2 (p 2 )(g(p 9 )) (2) 

and gluon-fusion process 

g{ Pg )W*(q) - gi(pi) ?2 (p2). (3) 

Hereafter, we shall assume that qi is light and q 2 is heavy, i.e. p\ = and p\ = m 2 . 
We also define Q 2 = —q 2 and the Bjorken x as 

* = w~ q ' (4) 

where P is momentum of the nucleon upon which the incoming lepton scatters, i.e. a 
proton for ep interaction, a proton or a neutron for uN processes. 

In the following, in order to appreciate quark mass effects, we shall be mainly inter- 
ested in values of Q 2 with respect to which m 2 is not negligible. Nevertheless, we shall 
present results for both small and large values of m 2 /Q 2 and comment on the transition 
between the massive case and the massless approximation. 
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At small values of Q 2 one should also take into account the effect of the target mass 
M, typically of the order of 1 GeV. This leads to the introduction of the Nachtmann 
variable rj, which generalizes the Bjorken x to the case of non-negligible target mass. It 
satisfies the following relation: 



I l_ /J_ NP_ 

r) ~ 2x + V 4x 2 + Q 2 ' U 

The coefficient function of the quark-scattering process (J2J) is usually expressed in terms 
of the variable 

Q 2 + m 2 

z = -7, > 6 

2pi • q 

with < z < 1. Here we also introduce the energy fraction £ of q\ with respect to the 
initial-state proton 1 and the quantity 

X = ~ y (7) 

where we have defined 

In the approximation of vanishing target and quark masses, \ coincides with the Bjorken 
x, i.e. x ~ > x - Moreover, the relation z = x/£ holds and x and \ are related via 



x = ^ f 9 ) 

1-M 2 X 2 X 2 /Q 2 ' 

The Bjorken x, which in the massless approximation varies between and 1, is now 
constrained in the range: 

0< ^ i-mW (10) 

The differential distributions da/dz of the parton- level processes in Eqs. (0) and (jSJ) 
are collinear divergent and their explicit expressions can be found in JH] • Such collinear 
singularities are usually regularized in dimensional regularization and subtracted in the 
MS factorization scheme. This procedure yields the MS coefficient functions [20J. 



Following the notation of [21], and referring to neutrino scattering on a fixed target, 
the inclusive cross section for viv) scattering in the CC regime can be expressed as a 
linear combination of three structure functions F\, F 2 and F 3 : 



d 2 o v{ ? ) GlME 



dxdy n(l + Q 2 /mlv) 2 



{y'xF, + [l - (l + ^) y] F 2 ± y (l - |) xF 3 } . (11) 



1 £ can be defined, e.g., in a frame where the nucleon and the quark are collinear |18j . Given 
P = (P(°),0,0,P( 3 )) and pi = (pi 0) ,0,0,pi 3) ) ) one has: £ - (p[° } + p[ 3) ) / (P (0) + P (3) ). 
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In the above equation, G F is the Fermi constant, E is the neutrino energy in the target 
rest frame and y the inelasticity variable, given by: 



y 



P-q 
P~k' 



:i2i 



with k being the momentum of the incoming lepton, the neutrino in the case of Eq. (jllj) . 
In Eq. (fTT|) . we have ignored structure functions F4 and F5 since their contribution to the 
differential cross section is proportional to powers of the lepton masses and is therefore 
suppressed [2T] • Analogous expressions to Eq. (fTTj) hold for positron (electron) scattering 
on a proton as well. 



Following 



one defines the 'theoretical' structure functions Ti which are 



related to the F^s of Eq. (fTTj) via the following relations: 



2x 

v 

2 
P 



T 



2- 



where we have defined: 



1 + 



(2Mx\ 



(13) 
(14) 

(15) 



(16) 



The introduction of the structure functions T is convenient as they can be straight- 
forwardly expressed as the convolution of parton distribution functions with MS coeffi- 
cient functions: 



T i(^c, Q 



1 d£ 

x T 



Ci(£,p 2 ,p 2 F ,\)qi 



|,/4j +C?{£,n 2 ,(i 2 F ,\)g lj,n 2 F 



(17) 



for i = 1,2,3. In Eq. (fTTj) . fi and fi F are the renormalization and factorization scales; 
Cf and Cf are the coefficient functions for quark-scattering (J2J) and gluon fusion (J3j) 
processes; qi(C, P>f) an d 9(C,Pf) are the initial-state light quark and gluon parton dis- 
tribution functions in (J2J) and (JHJ). We have chosen x as the argument of T since later 
we shall present results for the structure functions in terms of x. 

The coefficient functions can be evaluated perturbatively and, to 0(as), they read: 



C?(z,» 2 , t i F ,X) = 8(l-z) + 



2tt 



-Hf(z,ii 2 F ,X), 



C?(z : H 2 i11 2 f,\) 



2tt 



-Hf(z,fM 2 F ,X). 



(18) 
(19) 
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The expressions for the functions #1,2,3 with the collinear divergences subtracted in the 
MS factorization scheme can be found in pOj. For consistency, in Eq. (fTTj) one will have 
to use MS parton distribution functions as well. 

As we are inclusive with respect to the final-state heavy quark q2, the quark-scattering 
coefficient functions are free from mass logarithms ln(m 2 /Q 2 ). However, the functions 
Hf present terms which behave like ~ 1/(1 — z) + or ~ [hi(l ~ z )/{^ ~ z )]+ an d are 
therefore enhanced once the quark energy fraction z approaches one [201 • The z — > 1 
limit corresponds to soft-gluon emission. In Mellin space, such a behaviour corresponds 
to contributions ~ In iV or ~ IniV 2 in the limit N — > 00, where the Mellin transform of 
a function f(z) is defined as: 



N 



dzz N ~ 1 f{z). 



(20) 



The gluon-initiated coefficient functions Hf in Eq. (fTTIj) are however not enhanced in the 
limit z — > 1, since the splitting g — > qq is not soft divergent. 

We would like to make further comments on the coefficient functions at large z. 
Omitting terms that are not enhanced in the soft limit for any value of the mass ratio 
m 2 /Q 2 , one has: 

Hf(z^l^ 2 F ,X)=H soi \z,fx 2 F ,X), (21) 

where 



H soit (z,fi F ,\) 



+ 



2CV 

1 

4 



{• 


"ln(l-z)~ 




"ln(l - \z)~ 




1 - z 


+ 


1 - z 



1 - z 



1 - \z) 



+ 



In 



Q 2 + m 2 

/4 



(22) 



with C F = 4/3. 



One immediately sees that the behaviour of the coefficient function (j22j) at large z 
strongly depends on the value of the ratio m/Q. To illustrate this point, we consider 
the two extreme regimes: small-m/Q limit, i.e. A ~ 1, and large-m/Q limit, i.e. A 1. 



In the first case, setting A = 1 in (j22j) and performing its Mellin transform, one 
recovers the large- N limit of the iV-space MS massless coefficient function reported 
inHUl: 



Cjv (a* 1 A*F) Q 



A=l 



1 + 



a s (n 2 )C t 



71 



- In 2 iV + 
2 



IE + 



In 



/4 



IniV 



(23) 



where 7^ = 0.577 ... is the Euler constant. 
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In the second case, for A —>■ and z — *■ 1, one has Xz = A in (I22j) . As a result, the 
term ~ [(1 — z)/(l — \z) 2 ] + is regular in the soft limit. The following expression holds 
in moment space: 



c7V,/4,QW)Ia<<i = i + 



ln 2 iV + 



7T 



, .M 2 
2 7b + 1 - In — r 



In AT 



For later convenience, we have introduced the scale: 



(24) 



M' 



mf 1 + 



Q 2 



(25) 



Furthermore, in evaluating the Mellin transforms, we have made use of the following 
relation [B]: 



( l dzz N - 1 


~\n k (l-z)~ 




= 1 


10 







lnJV + 0(l) 



for = 0, 



l/2\n 2 N + -y E \nN + 0(l) for fc = 1. 



The different behaviour in the soft limit of the NLO MS coefficient function for small 
and large values of the ratio m/Q, which can be seen by comparing Eqs. (j23j) and ([24)1 . 
will lead to different expressions of the resummed coefficient function, which will depend 
on the value of m/Q. This will be discussed in detail in the next section. 



3 Soft-gluon resummation 

In this section we would like to discuss the resummation of soft-gluon radiation in the 
quark- initiated MS coefficient function of CC DIS. 

Soft resummation in DIS can be performed along the general lines of [8 : one evalu- 
ates the amplitude of the process qi(pi)W*(q) — > q2(P2)g(p g ) in the eikonal approxima- 
tion and exponentiates the result. 

However, we can think of a simpler solution to get the resummed coefficient function 
in massive CC DIS. We note that, in the eikonal approximation, the coefficient function 
for heavy quark production in charged-current DIS is related by crossing to the coefficient 
function of heavy-quark decay into a W and a light quark, for example top decay 
t(pt) W(q)b(pb)(g(p g )). In fact, in top decay the b quark can be treated in the 
massless approximation m& ~ 0, since m 2 -C m 2 and m 2 <C rriyy 

If in top quark decay one observes the b quark, i.e. computes the differential width 
with respect to the 6-quark energy fraction xj, \lb\ H7\. one recovers a situation which is 
analogous to CC DIS. We detect a light quark, namely the b in the top decay coefficient 
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function and gi in CC DIS, and are inclusive with respect to the heavy quark, the top 
quark in t-decay, the massive quark g 2 in process (J2J). 



The correspondence with the top-decay process implies that, for large N, the CC 
DIS large-m/Q coefficient function (J23|) coincides with the top decay coefficient function 
computed in Refs. after one performs the replacements: 



m 



m 2 w 



(26) 



We can therefore apply Eq. (|2fij) . along with p t —> P2 and p^ —> p% in the eikonal current 
of Ref. [T7j and repeat all the steps which have led to the resummation of the top decay 
MS coefficient function, which we do not report here for the sake of brevity, and resum 
soft contributions in CC DIS to next-to-leading logarithmic (NLL) accuracy. 

As in ^7], we write the resummed coefficient function as an integral over z and 
k 2 = (pi +Pg) 2 (l — z )i with /Mp < k 2 < Ai 2 (l — z) 2 , and obtain the following expression: 



In A 



N 



y N—l 



dz- 



- 1 



M\l-z? dk< 



A 



a s {k 2 ) +S a s [M 2 {l-z 



,(27) 



where M. has been introduced in Eq. (|25|). 



We point out that the upper limit of the integration over k 2 can be explicitly ob- 
tained, e.g., in the centre-of-mass frame, where the light-quark energy satisfies the rela- 
tion 

Q 2 + m 2 



2Ex 



(l-z)Q 2 + m 2 



28 



For finite m, in the soft limit z — > 1, 1E X ~ (Q 2 + m 2 )/m. From k 2 < AE 2 (1 — z) 2 , one 
gets the upper limit k 2 < M.{\ — z) 2 . 

The function A(cts) in Eq. (|27j) can be expanded as a series in as as: 



A {n) . 



(29) 



The first two coefficients are mandatory in order to resum our coefficient function to 
NLL accuracy jH]: 

A {1) = C F , (30) 

r2\ K 



A (2) 



1 



c f 



C A 



67 
18 



7T" 



9 



(31) 



where Cf = 4/3, Ca = 3 and n/ is the number of quark flavours. A numerical estimate 
for is also known [2H] (see also [TT]). 
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The function S(cts) is characteristic of processes where a heavy quark is involved 
and expresses soft radiation which is not collinear enhanced. 



The following expansion holds: 



n=l 



where only the first coefficient is needed to NLL level. It reads: 



(32) 



(33) 



The integral over z in Eq. (J27|) can be performed, up to NLL accuracy, by making the 
following replacement [Hj: 



z N - 1 - 1 



-e i 



"7-B 



N 



(34) 



where is the Heaviside step function. This leads to writing the following result for 
the function A N : 



A N (fi 2 , n%, m 2 , Q 2 ) = exp \\nNg^(£) + g {2 \£, /i 2 , 



with 



I = b a s (n 2 )lnN 
and the functions and gw given by 

A {i) 



s ( W,/4) 



[2£ + (l-2£) ln(l-2£)] 



In 



M 2 

f4 



2 7i 



ln(l - 2£) 



+ 



4£ + 21n(l-2f)+ln 2 (l-2f) 



In Eqs. 



2nb i 
AV_ 
2nb 

-^[2£ + ln(l-2£)]f^ + A«ln4 
2tt6 \7r6o /4, 

+ —-In (l-2£). 
2nbn 



bo and b\ are the first two coefficients of the QCD ft function: 

33 — ri/- , 153 — 19nf 
b o = 77; h = 



(35) 
(36) 

(37) 



12vr 



24tt 2 



(3* 



(39) 



In Eq. (135)1 the term In Ng^'(£) accounts for the resummation of leading logarithms 
(LL) ag\n n+1 N in the Sudakov exponent, while function g^(£ , /ii 2 , [i 2 F ) resums NLL 
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terms In™ N. One can also check that the 0(as) expansion of Eg. ()85jl reproduces the 
large- N coefficient function in Eq. (J2IJ). 

Furthermore, we follow ^HE] and also include in our final Sudakov-resummed 
coefficient function terms which, in the iV-space coefficient function, are constant with 
respect to iV and which we denote with K^/j^, m 2 , Q 2 ): 



HUfi 2 ,fi 2 F ,m 2 ,Q 2 ) 



1 + 



2tt 



■/Q(/4,m 2 ,Q 2 ) 



x exp 

The explicit expression for the constant terms reads 

' Q 2 + m 2 \ 



\nNgU{e)+g®{e,v,LiF) 



(40) 



3 
2 " 



hi 



[l 2 F 



+ ln(l - A) 2 lE 



3A-2 X 



2A 



+ A-2Li 2 - 



Q 2 



+ 2( 7e -1)( 7e + 2), 



where A was introduced in Eq. 
given by: 



and we have defined the quantities Ai 



A 2 = -(1-A)ln(l-A). 



(41) 
which are 

(42) 
(43) 



We now match the resummed coefficient function to the exact first-order result: we 
add the resummed result to the exact coefficient function and, in order to avoid double 
counting, we subtract what they have in common. Our final result for the resummed 
coefficient function reads: 

tf^(/i 2 ,/4,m 2 ,Q 2 ) = ^(/i 2 ,/4,m 2 ,Q 2 ) 

- H Nii (n 2 ,n 2 F ,m 2 ,Q 2 )^ a 
+ Cjy/i 2 ,/4,m 2 ,Q 2 ) S (44) 
where C q Ni is the Mellin transform of the exact 0(as) coefficient functions Cf in Eq. (JTBJ). 

Before closing this section, we would like to compare our calculation with the results 
on soft resummation in DIS processes when the final-state quark g 2 is light |KJEIIJE]- 
In fact, while for top decay, up to width effects, the value of mf/m^ is fixed and rather 
large, m 2 /Q 2 can vary widely in CC DIS according to the value of Q 2 in the considered 
process. This implies that, if m 2 <C Q 2 , the heavy quark behaves in a massless fashion 
and one may have to use a different result for the resummed coefficient function. 

As we have already remarked in Section 2, one cannot just recover the massless limit 
from the massive result: the limits m/Q — > or A — » 1 of Eqs. (|2T|) and (J4T!|) appear 
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in fact to be divergent. The actual result is instead finite: one cannot naively apply 
those equations for very small values of m/Q, but rather needs to perform a different 
calculation in this limit. 

In fact, if the final-state q2 is massless, the NLL resummed coefficient function re- 
ceives an additional contribution due to purely collinear radiation. As shown in [B] , one 
can account for such a contribution by solving a modified evolution equation for the 
parton distribution associated with the final state. 

Therefore, in the massless approximation, one will have to add an additional term 
to S(atg) to account for collinear-enhanced radiation, which leads to the introduction 
in the Sudakov exponent of a function that is usually denoted by l/2B[a s (Q 2 (l — z))). 
The argument of function B(as) is indeed the virtuality of the final-state jet, which, for 
a massless qi in process (J2J), will vanish in soft limit like (j>2 + p g ) 2 — (1 — z)Q 2 . 

On the contrary, if the final-state quark is massive, the virtuality of the final state in 
the large- z limit is given by (p2+p g ) 2 = (l—z)(Q 2 +m 2 )+m 2 and is never small, provided 
m/Q is sufficiently large. We do not have therefore additional collinear enhancement for 
heavy quark production, as we do not have any ~ B(as) contribution in the resummed 
top decay coefficient function of Ref. [T7j . 

Since in the following section we shall present results in the limit m/Q — > as well, 
we explicitly write the resummed coefficient function in the massless approximation . 
As for the large-m/Q case, we write it as an integral over z and k 2 = (pi + p g ) 2 (l — z), 
with/4< k 2 <Q 2 (l-z): 




2(1-*) ^2 

lnAjvU/Q-o = I dz~— — \ I o -j^-A 



a s (k 2 )} + l -B [a s (Q 2 (l 



(45) 

The upper limit for the integration variable k 2 can be obtained using Eq. (j2HI) : for 
m = and z — *■ 1, one has 2E\ ~ Q/yl — z. From k 2 < 4Ef(l — z) 2 , the relation 
k 2 < Q 2 (l - z) follows. 

We can expand function B(as) as a series in a^: 

B(as) = f(^-) n B^ (46) 



7T 



n=l 

and, to NLL level, keep only the first term of the expansion jH]: 

= ~\C F . (47) 
Explicit results for Eq. (|43j) to NLL level can be found in Ref. [TT], where functions 
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analogous to our and g^ 1 are reported 2 . One can also check that the 0(as) 
expansion of Eq. (|45|) yields the large- N limit of the NLO coefficient function presented 
in Eq. 

As in Eq. (|4()jl . we also include in the final Sudakov-resummed massless coefficient 
function terms which are constant with respect to N. Once again, we observe that 
such constant terms cannot be obtained just taking the m/Q — > limit of the massive 
case, i.e. Eq. (j41j) . but we shall have to explicitly extract them from the massless NLO 
coefficient function. We obtain: 

^G4,m 2 ,Q 2 )| m/( ^ = g_2 7JS )lnP^ + 7 | + \lE - y - \, (48) 

for i = 1, 2, 3. As in Eq. we match the resummed result (|4"5|) to the exact coefficient 
function ()18|l. with functions iff now evaluated in the limit m/Q — ► 0, i.e. A — > 1. 

More details on the derivation of Eq. ()27|) as well as on the comparison with the 
massless result fl43j) can be found in |27| . 



4 Phenomenological results 

We would like to present results on proton structure functions for charm quark pro- 
duction in CC DIS and investigate the effects of soft-gluon resummation. As discussed 
in Section 2, structure functions are given by convolutions of MS coefficient functions 
with parton distribution functions. Since the resummed coefficient function is given in 
iV-space, in principle one may like to use parton distribution functions in Mellin space 
as well, in order to get the resummed structure function in moment space and finally in- 
vert it numerically from N- to x-space. However, all modern sets of parton distribution 
functions [21112111201 are given numerically in the form of grids in the (x, Q 2 ) space. 

In order to overcome this problem, we follow the method proposed in Ref. j2H| in 
the context of joint resummation, which allows one to use x-space parton distributions, 
even when performing resummed calculations in Mellin space. 

Such a method consists of rewriting the integral of the inverse Mellin transform of 
the resummed structure functions as follows: 

fr(x, Q 2 ) = ^- J dN X - N qN^ 2 F )H^(fi 2 , m 2 , Q 2 ) 

2 We point out that Ref. fTJ contains results to next-to-next-to-leading logarithmic (NNLL) accuracy 
for soft resummation in massless DIS. However, for the sake of comparison with the massive result i|27[l . 
we shall use in the following the results of to NLL. 
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1 f J AT ~N Wat 1\2 ! 2\^N%V 2 ^ 2 F-> m2 ->Q 2 

2^L dNX [C^- 1 )^)] 1^-2 



where qNifJ^) is the parton distribution of the initial-state quark qi, as if it were available 
in N space, and T^v is the integration contour in the complex plane, chosen according 
to the Minimal Prescription [26J. 

One can integrate out the term with the parton distribution function by observing 
that the following relation holds: 



JL ^ dm -« (N _ d w ,) = | {4 [&«, ^,)] } = (50) 



Due to the above relation, one does not need anymore the iV-space parton distributions, 
but one should just differentiate the x-space ones, which numerically can be done. 

As one has the analytical expression for in iV-space (J44j) . the following inverse 
Mellin transform is straightforward: 



w(e,^) = T L r/ dNc N ( 5i ) 



2ni Jr N 5 (N - l) 5 

which, thanks to the suppressing factor 1/(N — l) 2 in the integrand, turns out to be 
smooth for £ — > 1. 

The resummed structure function will be finally expressed as the following convolu- 
tion: 



(52) 



After having clarified how to deal with the parton distribution functions, we are 
able to present results for soft-resummed structure functions. We shall consider charm 
quark production, i.e. q 2 = c in Eqs. (JTJ) and (J2J), since processes with charm quarks in 
the final state play a role for structure function measurements at NuTeV and HERA 
experiments. 

As for the choice of the parton distribution set, in principle, once data on heavy 
quark production in CC DIS is available, one should use the NLL coefficient functions 
when performing the parton distribution global fits and get NLL parton densities as 
well. For the time being, however, we can just use one of the most-updated NLO sets 
and convolute it with the fixed-order or resummed coefficient function. 

We shall present results based on the new generation of CTEQ NLO MS parton dis- 
tribution functions |2H], the so-called CTEQ6M set, but similar results can be obtained 
using, e.g., the MRST |2H] or the GRV [SO] sets. 
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The elementary scattering processes which yield the production of charm quarks in 
CC DIS are dW* — > c and sW* — > c. For e + (e~)p scattering at HERA, our parton 
distribution function qi((,,Q 2 ) in Eq. (JT7J) will be as follows: 

<?i(£> Q 2 ) I hera = \Vcd\ 2 d(Z,Q 2 ) + \V^s{^Q 2 ), (53) 

where V c d and V cs are the relevant Cabibbo-Kobayashi-Maskawa matrix elements. For 
neutrino scattering on an isoscalar target, Eq. (j53j) will have to be modified in order to 
account for the possibility of an interaction with a neutron as well: 

gi(e,Q 2 )|NuTev = \ Vcd \ ^ Q2)+ 2 U ^ Q2) + \V cs \ 2 s^,Q 2 ), (54) 

where d(^,Q 2 ), u(^,Q 2 ) and s(£,Q 2 ) are still the proton parton distribution functions 
and we have applied isospin symmetry, i.e. u p = d n and s p = s n . 

In order to be consistent with the use of the CTEQ parton distribution functions, we 
shall use for the QCD scale in the MS renormalization scheme the values A 4 = 326 MeV 
and A5 = 226 MeV, for four and five active flavours respectively. This corresponds to 
as{mz) — 0.118. The charm and bottom quark masses have been set to m c = 1.3 GeV 
and mi, = 4.5 GeV, as was done in [2*%] . 

Since mass effects are important for large values of the ratio m 2 /Q 2 , we would like 
to have Q as small as possible so as to be able to apply our massive resummation. The 
NuTeV experiment is indeed able to measure structure functions at small Q 2 values 
[TJI2], as it can detect the final-state muon produced in v^N — > fiX processes. In our 
phenomenological analysis we shall consider charm production at Q 2 = 2 GeV 2 and 
Q 2 = 5 GeV 2 , which are values reachable at NuTeV such that m 2 /Q 2 is relatively large. 

The detection of CC events at HERA is more problematic, due to backgrounds and 
to the presence of a neutrino, instead of a charged lepton, in the final state. The current 
high-luminosity HERA II experiments may detect heavy quarks in CC events, but the 
Q 2 values for such events are still supposed to be much larger than the charm quark 
mass squared, typically Q 2 > 100 GeV 2 [31] . We shall nonetheless present results for 
charm quark production at HERA, but in this case we shall have to use the massless 
result for the resummed MS coefficient function, i.e. Eq. (J4*5|) and the formulas reported 
in [IT]. We shall consider the typical HERA values Q 2 = 300 GeV 2 and 1000 GeV 2 . 

We present results for the structure function F£(x, Q 2 ), but we can already anticipate 
that the effect of the resummation is approximately the same on all three structure 
functions Ff, F| and F| and on the single-differential cross section da/dx, obtained 
from Eq.([lip. after integrating over y. 

For the sake of comparison with the experiments, we shall plot the structure function 
_F 2 C in terms of the Bjorken x, which is the measured quantity. For most of our plots, 
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Figure 1: Results for F^ix^Q 2 ) for charm quark production in neutrino scattering at Q 2 = 
2 GeV 2 with (solid) and without (dashed) soft resummation in the coefficient function. We 
have set \xp = /i = Q. 

we shall set renormalization and factorization scales equal to Q, i.e. fip = fi = Q. 
Afterwards, we shall also investigate the dependence of our results on the choice of such 
scales. 
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Figure 2: As in Fig.HJ but for Q 2 = 5 GeV 2 . 
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Figure 3: Results for F| (x, Q 2 ) for positron-proton scattering at Q 2 = 300 GeV 2 , with (solid) 
and without (dashes) soft resummation in the coefficient function. We have set \ip = [i = Q. 
In the inset figure, we show the same plots at large x and on a linear scale. 

In Figs. Qand El we show F£(x, Q 2 ) in neutrino scattering at Q 2 = 2 GeV 2 and 5 GeV 2 
and compare the results obtained using fixed-order and soft-resummed MS coefficient 
functions. We have set the target mass to M — 1 GeV, which is a characteristic nucleon 
mass. We observe a relevant effect of the implementation of soft resummation: the two 
predictions agree up to x ~ 0.2 — 0.3, afterwards one can see an enhancement of the 
structure function due to the resummation. At very large x, resummed effects are indeed 
remarkable: for Q 2 = 2 GeV 2 one has an enhancement of a factor of 2 at x = 0.5 and 
a factor of 5 at x = 0.6. For Q 2 = 5 GeV 2 one gets a factor of 2 at x = 0.7 and 5 at 
x = 0.8. 

In Figs. El and 0] we present results for F£(x, Q 2 ), but for charm production at 
HERA, in particular for positron-proton scattering at Q 2 = 300 and 1000 GeV 2 . In 
this case, since m c /Q C 1, we use the massless result f|4*5]) for the resummed coefficient 
function. We see that the impact of the resummation is smaller than in the case of 
low Q 2 values. This is a reasonable result: in fact, leading and next-to- leading loga- 
rithms in the Sudakov exponent are weighted by powers of «s(yU 2 ). As, for example, 
as (2 GeV 2 ) ~ 3 as (300 GeV 2 ), resummed effects are clearly more important when Q 2 is 
small. Moreover, the larger Q 2 is, the larger the values of x are at which one is sensitive 
to Sudakov effects. 

We note in Figs. © and (HJ) that, for x > 0.6, fixed-order and resummed predictions 
start to be distinguishable. We estimate the overall impact of large- x resummation on 
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Figure 4: As in Fig.EJ but for Q 2 = 1000 GeV 2 . 



F 2 c (x, Q 2 ) at Q 2 = 300 GeV 2 and Q 2 = 1000 GeV 2 to be between 10% and 20%. 

In Figs. El and El we investigate the dependence on factorization and renormalization 
scales. We still keep /x = /ip, but we allow such scales to assume the values Q 2 , 2Q 2 and 
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Figure 5: Dependence of F% on the factorization scale for neutrino scattering at Q 2 = 2 GeV 2 . 
Solid lines include soft resummation in the coefficient function, dashed lines are fixed-order 
predictions. 
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Figure 6: As in Fig. [51 but for e + p scattering at Q 2 = 300 GeV 2 . In the inset figure, we plot 
the same curves at large x and on a linear scale. 

AQ 2 . We consider Q 2 values of 2 and 300 GeV 2 for the experimental environments of 
NuTeV and HERA respectively. We see that the curves which implement soft resum- 
mation in the coefficient function show a weaker dependence on the chosen value of the 
factorization/renormalization scales. In Fig. |S]one can see that one still has a visible ef- 
fect of the value of such scales, but the overall dependence on /x F and \i of the resummed 
prediction is smaller than for the fixed-order. The plots at large Q 2 exhibit, in general, a 
weak dependence on fip and fi even at NLO, as shown in Fig. in fact, the dependence 
on the factorization and renormalization scales is logarithmic, and hence smaller once \i 
and fip vary around large values of Q. Nevertheless, while the NLO structure function 
still presents a residual dependence on the scales, the resummed predictions obtained 
with three different values of //p are basically indistinguishable. A smaller dependence 
on such scales implies a reduction of the theoretical uncertainty of the prediction and is 
therefore a remarkable effect of the implementation of soft gluon resummation. 

We finally would like to compare the impact that mass effects and soft-gluon resum- 
mation have on the charm structure functions. To achieve this goal, we plot in Fig. the 
theoretical structure function (x, Q 2 ) defined in Eq. ([T7]) as a function of the variable 
X (see Eq. 0), for neutrino scattering at Q 2 = 2 GeV 2 . We compare fixed-order massive 
(dashed line), fixed-order massless (dots) and massive resummed (solid) predictions. We 
observe that the two fixed-order calculations yield different predictions throughout the 
full x range, which is a consequence of the implementation of mass effects; however, 
at large Xi where one starts to be sensitive to the resummation, the impact of soft 
resummation is competitive with the one of mass contributions. 
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Figure 7: Comparison of massive (dashes) and massless (dots) fixed-order calculations with 
the resummed massive result (solid) for Q 2 = 2 GeV 2 and fi = (j,f = Q- Plotted is the 
theoretical structure function ^(.XiQ 2 ) defined in Eq. (|T7|) . 

Before closing this section, we would like to point out that, although we have im- 
proved our perturbative prediction implementing soft-gluon resummation in the coeffi- 
cient function, non-perturbative corrections are still missing. Non-perturbative effects 
are important especially at small values of Q 2 and large x. In fact, a more accurate 
investigation of the very-large x limit of the structure functions shows that they exhibit 
an oscillatory behaviour once x gets closer to the maximum value which is kinematically 
allowed. 



5 Conclusions 

We have studied inclusive heavy quark production in charged-current Deep Inelastic 
Scattering. We have reviewed the results for the NLO cross section and structure func- 
tions and observed that the MS quark-initiated coefficient function contain terms which 
become arbitrarily large once the initial-state light-quark energy fraction approaches 
unity, which corresponds to soft-gluon radiation. We have resummed soft-gluon contri- 
butions in the MS coefficient function to next-to-leading logarithmic accuracy. We have 
compared the results with the resummed coefficient function for 6-quark production in 
top decay and for CC DIS processes, but with a light quark in the final state. 

We have presented results on the structure function F£ (x, Q 2 ) for charm quark pro- 
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duction at NuTeV and HERA, using the massive or the massless coefficient function 
according to the value of the ratio m c /Q. We have compared predictions obtained con- 
sidering fixed-order and resummed coefficient functions and have found a remarkable 
effect due to soft resummation. The structure functions exhibit an enhancement at 
large x, which is visible especially when Q is comparable with the heavy quark mass. 
The results including soft resummation exhibit very little dependence on the choice of 
factorization and renormalization scales, which is a reduction of the theoretical uncer- 
tainty Furthermore, we have compared massive and massless predictions at small Q 2 
and have found that, at large x, resummation effects are of a magnitude similar to mass 
corrections. 

It will be clearly very interesting to compare our predictions with experimental data. 
New data on NuTeV structure functions will be soon available j22] and we plan to 
investigate whether the structure functions yielded by our calculation are able to fit the 
data. In particular, we wish to evaluate the role played by soft-gluon resummation in 
performing such fits. The preliminary analyses in Ref. pQ show in fact that the NuTeV 
data lies above the NLO structure functions based on Ref. [33] for x > 0.55. The effect 
due to the resummation is indeed an enhancement of structure functions at large x, 
hence it may throw some light on the NuTeV studies. However, for such a comparison 
to be reliable, we shall have to implement nuclear-correction effects into our structure 
functions, to account for the neutrino scattering on an iron target. A similar comparison 
can be made with the foreseen CC data from HERA II. Moreover, we believe that once 
such a data becomes available, our calculation can also be used to perform updated NLL 
global fits of the parton distribution functions. 
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